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We discuss the O(a) electroweak radiative corrections to W boson production at the Tevatron and their effect on the W boson 
mass extracted by experiment. The results of a new calculation of the O(a) corrections are presented and compared with those of a 
previous calculation. We also briefly discuss the O(a) corrections to Z boson production at the Tevatron and two- photon radiation 
in W and Z events. 



1 Introduction 

The Standard Model of electroweak interactions (SM)| 
so far has met all experimental challenges and is now 
tested at the 0.1% levelEl. However, there is little direct 
experimental information on the mechanism which gen- 
erates the masses of the weak gauge bosons. In the SM, 
spontaneous symmetry breaking is responsible for mass 
generation. The existence of a Higgs boson is a direct 
consequence of this mechanism. At present the nega- 
tive result of direct searches performed at LEP2 imposes 
a lower bound of M H > 89.8 GeV0 on the Higgs bo- 
son mass. Indirect information on the mass of the Higgs 
boson can be extracted from the Mh dependence of ra- 
diative corrections to the W boson mass, Mw, and the 
effective weak mixing angle, sin 2 O^j. Assuming the SM 
to be valid, a global x 2_n t to all available electroweak 
precision data yields a 95% confidence level upper limit 
on M H of 280 GeVH. 

The current estimate of Mh strongly depends on 
the world average for the weak mixing angle, sin 2 6^¥l = 
0.23155±0.00019i. It results from a combination of LEP 
and SLC data which currently are not in good agree- 
ment □. Furthermore, sin 2 O^ff 1S Quite sensitive to the 
hadronic contribution to a(M§), Aai ia d(M§). The ac- 
curacy of Aahad(Af|) has been the subject of a num- 



ber of publications during the last four years E 
ror estimates range between <5(Aahad(-W§)) = 0.0007 
and 5(Aa ha d(M|)) = 0.00016 B. A smaller error for 
Aahad(M§) implies that sin 2 6^Fl receives more weight 
in the Mh fit, i.e. the discrepancy between the LEP and 
SLC data becomes a limiting factor in the estimate of 
the Higgs boson mass from electroweak data. 

A more precise measurement of Mw is, therefore, 
very important in order to extract more accurate infor- 
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mation on Mh from electroweak data. Furthermore, in 
contrast to sin 2 , the W mass depends only mildly on 
Aahad(Afi)i- A more precise measurement of Mw thus 
automatically reduces the sensitivity of the extracted 
Higgs boson mass to Aai, a d(M|) 1 _ Currently, the W bo- 
son mass is known to ±0.06 GeVtl from direct measure- 
ments. A significant improvement in the W mass uncer- 
tainty is expected in the near future from measurements 
at LEP20 and the Tevatron El. The ultimate precision 
expected for Mw from the. combined LEP2 experiments 
is approximately 40 MeVQ. At the Tevatron, integrated 
luminosities of order 2 fb _1 are foreseen for Run II, and 
one expects to measure, the W mass with a precision of 
approximately 40 Me Veil per experiment and decay chan- 
nel. 

In order to measure the W boson mass with high 
precision in a hadron collider environment, it is neces- 
sary to fully understand and control higher order QCD 
and electroweak (EW) corrections to W production. The 
determination of the W mass in a hadron collider envi- 
ronment requires a simultaneous precision measurement 
of the Z boson mass, Mz, and width, Tz- When com- 
pared to the value measured at LEP, the two quantities 
help to accurately determine the energy scale and resolu- 
tion of the electromagnetic calorimeter, and to constrain 
the muon momentum resolution!! It is therefore also nec- 
essary to understand the higher order EW corrections to 
Z boson production in hadronic collisions. 

Recently, new and more accurate calculations of the 
O(a) EW corrections—to 13 and Z boson production 
in hadronic collisions Hj became available. In a previous 
calculation, only thiLfinal state photonic corrections were 
correctly included lid. The sum of the soft and virtual 
parts was estimated from the inclusive 0(a 2 ) W — > lv{"j) 
and Z — » £ + £~(-f) (I = e, fi) width and the hard photon 
bremsstrahlung contribution. Initial state, interference, 
and weak contributions to the 0(a) corrections were ig- 
nored altogether. The unknown part of the 0(a) EW 
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corrections in Ref. [jnj, combined with effects of multi- 
ple photon emission, have been estimated to contribute 
a systematic uncertainty of 8Mw = 15 — 20 MeV to the 
measurement of the W massB. 

In Section 2, we briefly describe the technical details 
of the calculation of the 0(a) corrections to W boson 
production presented in Ref. [||, and compare the re- 
sults with those of Ref. [[fl]] . In Section 3 we summarize 
the calculation of the 0(a) QED corrections to Z boson 
production reported in Ref. ]l0| , and in-Section 4 some 
preliminary results of a new calculation 113 of two-photon 
radiation in W and Z production in hadronic collisions 
are presented. 



a) pp^e*v{j) 
Vs = 1.8 TeV 
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2 Electroweak Corrections to W Boson Produc- 
tion at the Tevatron 

The calculation of the 0(a) corrections to pp — > W — > 
ivu is based on the full set of 0(a 3 ) Feynman diagrams, 
and includes both initial and final state radiative correc- 
tions, as well as the contributions from their interference. 
Final state charged lepton mass effects arc included in the 
following approximation. The lepton mass regularizes 
the collinear singularity associated with final state pho- 
ton radiation. The associated mass singular logarithms 
of the form hi(s/mf), where s is the squared parton cen- 
ter of mass energy and me is the charged lepton mass, 
are included in our calculation, but the very small terms 
of 0(mf/s) are neglected. 

To perform the calculation, a Monte Carlo method 
for next-to-leading-order (NLO) calculations similar to 
that described in Ref. E3] was used. With the Monte 
Carlo method, it is easy to calculate a variety of observ- 
ables simultaneously and to simulate detector response. 
Calculating the EW radiative corrections to W boson 
production, the problem arises how an unstable charged 
gauge boson can be treated consistently in the frame- 
work of perturbation theory. This problem has been 
studied in Ref. [ |l4"[ with particular emphasis on finding 
a gauge invariant decomposition of the EW 0(a) correc- 
tions into a QED-like and a modified weak part. In W 
production, the Feynman diagrams which involve a vir- 
tual photon do not represent a gauge invariant subset. In 
Ref. jl| it was demonstrated how gauge invariant con- 
tributions that contain the infrared (IR) singular terms 
can be extracted from the virtual photonic corrections. 
These contributions can be combined with the also IR- 
singular real photon corrections in the soft photon region 
to form IR-finite gauge invariant QED-like contributions 
corresponding to initial state, final state and interference 
corrections. The IR finite remainder of the virtual pho- 
tonic corrections and the pure weak one-loop corrections 
can be combined to separately gauge invariant modified 
weak contributions to the W boson production and decay 



Figure 1: Ratio of the O(o 3 ) and lowest order cross sections as a 
function of the transverse mass for a) pp — > e + !/(7) and b) pp — > 
fl^u('y) at y/s = 1.8 TeV for various individual contributions. The 
upper (lower) solid lines show the result for the QED-like initial 
(final) state corrections. The upper (lower) dotted lines give the 
cross section ratios if both the QED-like and modified weak initial 
(final) state corrections are included. The dashed lines display the 
result if only the initial - final state interference contributions are 
included. 

processes. 

The collinear singularities associated with initial 
state photon radiation can be removed by universal 
collinear counter terms generated by "renormaUzing" the 
parton distribution functions (PDF's) Ej, in complete 
analogy to gluon emission in QCD. In addition to the 
collinear counterterms, finite terms can be absorbed into 
the PDF's, introducing a QED factorization scheme de- 
pendence. We have carried out our calculation in the 
QED DIS and QED MS scheme. In order to treat the 
0(a) initial state QED-like corrections to W production 
in hadronic collisions in a consistent way, QED correc- 
tions should be incorporated in the global fitting of the 
PDF's using the same factorization scheme which has 
been employed to calculate the cross section. Current 
fits to the PDF's do not include QED corrections. A 
study of the effect of QED corrections on the_evolution 
of the parton distribution functions indicateslia that the 
modification of the PDF's is small. The missing QED 
corrections to the PDF introduce an uncertainty which, 
however, is likely to be smaller than the present uncer- 
tainties on the parton distribution functions. 

Since hadron collider detectors cannot directly de- 
tect the neutrinos produced in the leptonic W boson 
decays, W — > Iv, and cannot measure the longitudinal 
component of the recoil momentum, there is insufficient 
information to reconstruct the invariant mass of the W 
boson. Instead, the transverse mass (Mr) distribution 
of the Iv system is used to extract M\y. The various 
individual contributions to the EW O(o) corrections of 
the Mr distribution are shown in Fig. H. To compute the 
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cross section, we have_u.sed here the MRSA set of parton 
distribution functionsEJl The detector acceptance is sim- 
ulated by imposing the following transverse momentum 
(j>t) and pseudo-rapidity (77) cuts: 



p T (£)>25GeV, \i](t)\ < 1.2, 
i> T > 25 GeV. 



e, n, 



(1) 

(2) 



These cuts approximately model the acceptance cuts 
used by the CDF and D0 collaborations in their W mass 
analyses. Uncertainties in the energy and momentum 
measurements of the charged leptons in the detector are 
simulated in the calculation by Gaussian smearing of the 
particle four-momentum vector using—the specifications 
for the upgraded Run II D0 detector EH. 

The initial state QED-like contribution uniformly 
increases the cross section by about 1% for electron 
(Fig. [l|a,) and muon (Fig. [l]b) final states. It is largely 
canceled by the modified weak initial state contribution. 
The interference contribution is very small. It decreases 
the cross section by about 0.01% for transverse masses 
below My/, and by up to 0.5% for Mt > My/. The fi- 
nal state QED-like contribution significantly changes the 
shape of the transverse mass distribution and reaches 
its maximum effect in the region of the Jacobian peak, 
Mt ~ M\v- Since the final state QED-like contribution 
is proportional to ln(s/m^), its size for muons is con- 
siderably smaller than for electrons. As for the initial 
state, the modified weak final state contribution reduces 
the cross section by about 1%, and has no effect on the 
shape of the transverse mass distribution. 

In Fig. |l|, we have not taken into account realistic 
lepton identification requirements. When these require- 
ments are included, the mass singular logarithmic terms 
are eliminated in the electron case because the electron 
and photon momentum four vectors are combined for 
small opening, angles where it is difficult to resolve the 
two particlesLi. This significantly reduces the size of the 
EW corrections. On the other hand, in order to exper- 
imentally identify muons, the energy of the photon is 
required to be smaller than a critical value if the fj, — 7 
separation is small, and mass singular terms survive. Re- 
moving energetic photons thus enhances the effect of the 
O(a) corrections, and the effect of the EW corrections 
in the muon case is larger than in the electron case once 
lepton identification requirements are included. 

As we have seen, final state bremsstrahlung has a 
non-negligible effect on the shape of the Mt distribution 
in the Jacobian peak region. It is well known that EW 
corrections must be included when the W boson mass is 
extracted from data, otherwise the mass is shifted to a 
lower value. In the approximate treatment of the elec- 
troweak corrections used so far by the Tevatron exper- 
iments, only final state QED corrections are taken into 
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Figure 2: Ratio of the Mt distributions obtained with the QED- 
like final state correction p art of Rcf. H] to the one obtained using 
the approximation of Ref. for pp — > t^~v{-y) at yfs = 1.8 TeV. 



account; initial state, interference, and weak correction 
terms are ignored. Furthermore, the effect of the final 
state soft and virtual photonic corrections is estimated 
from the inclusive 0(a 2 ) W — > lv(ff) .width and the hard 
photon bremsstrahlung contribution Hi. When detector 
effects are included, the approximate calculation leads 
to a shift of about —50 MeV in the electron case, and 
approximately —160 MeV in the muon caseQ. 

Initial state and interference contributions do not 
change the shape of the Mt distribution significantly 
(see Fig. [lj) and therefore have little effect on the ex- 
tracted mass. However, correctly incorporating the final 
state virtual and soft photonic corrections results in a 
non-negligible modification of the shape of the transverse 
mass distribution for Mt > My/. This is demonstrated 
in Fig. [2[ which shows the ratio of the Mt distribution 
obtained with the QED-like final state correction part of 
our calculation to the one obtained using the approxima- 
tion of Ref. 0. 

The difference in the line shape of the Mt distribu- 
tion between the 0(a 3 ) calculation of Ref. || and the 
approximation used so far occurs in a region which is 
important for both the determination of the W mass, 
and the direct measurement of the W width. The pre- 
cision which can be achieved in a measurement of My/ 
using the transverse mass distribution strongly depends 
on how steeply the Mt distribution falls in the region 
Mt ~ My/. Any change in the theoretical prediction 
of the line shape thus directly influences the W mass 
measurement. From a maximum likelihood analysis the 
shift in the measured W mass due to the correct treat- 
ment of the final state virtual and soft photonic correc- 
tions is found to be AM W w O(10 MeV). This shift 
is much smaller than the present, uncertainty for My/ 
from hadron collider experiments a. However, for future 
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Figure 3: Ratio of the 0(a 3 ) and lowest order differential cross 
sections as a function of the di-lepton invariant mass for pp — > 
£+£-(1) at v / s = 1.8 TeV. 

precision experiments, a difference of 0(10 MeV) in the 
extracted value of M\y can no longer be ignored, and the 
complete 0(a 3 ) calculation should be used. 

3 Electroweak Corrections to Z Boson Produc- 
tion at the Tevatron 



lepton identification requirements are taken into account 
in Fig. |^. When these are included, the maximum en- 
hancement is reduced to a factor 1.6 for both electrons 
and muons. As for the W, initial state corrections are 
uniform and small. Final state radiative corrections dom- 
inate over the entire di-lepton invariant mass range. 

From Fig. || it is clear that final state bremsstrahlung 
severely distorts the Breit-Wigner shape of the Z reso- 
nance curve. As a result, QED corrections must be in- 
cluded when the Z boson mass is extracted from data, 
otherwise the mass extracted is shifted to a lower value. 
As in the W case, in the approximate treatment of the 
QED corrections to Z boson production used so far by 
the Tevatron experiments, only final state corrections 
are taken into account, and the effects of soft and vir- 
tual corrections are estimated from the inclusive 0(a 2 ) 
Z — > £ + £~ (7j)_width and the hard photon bremsstrahlung 
contribution til. When detector effects are taken into ac- 
count, the approximate calculation leads to a shift of the 
Z mass of about —150 MeV in the electron case, and ap- 
proximately —300 MeV in the muon caseQ. The Z boson 
mass extracted from our 0(a 3 ) £ + £~ invariant mass dis- 
tribution is found to be about 10 MeV smaller than that 
obtained using the approximate calculation of Ref . jll[ . 
This translates into an additional shift of several MeV in 
Mw through the dependence of the energy scale and the 
momentum resolution on the Z boson mass measured. 



The calculation of the 0(a) corrections to Z boson pro- 
ductionlilJ employs the same Monte Carlo method which 
was used in the W case. The collinear singularities orig- 
inating from initial state photon radiation are again re- 
moved by counter terms generated by renormalizing the 
PDF's. However, in contrast to W production, the Feyn- 
man diagrams contributing to the O(a) QED corrections 
can be separated into gauge invariant subsets correspond- 
ing to initial and final state corrections. Furthermore, the 
purely weak corrections form a separately gauge invari- 
ant set of diagrams. The weak corrections are expected 
to be very small and are therefore ignored in our calcu- 
lation. Both Z and photon exchange diagrams with all 
7 — Z interference effects are incorporated. 

In Fig. H we display the ratio of the 0(a 3 ) and the 
Born cross section as a function of the £ + £~ invariant 
mass in pp -> Z -> £+£~. For 40 GeV < m(l+l') < 
110 GeV, the cross section ratio is seen to vary rapidly. 
Below the Z peak, QED corrections enhance the cross 
section by up to a factor 2.7 (1.9) for electrons (muons). 
The maximum enhancement of the cross section occurs 
at m[i + £~) « 75 GeV. At the Z peak, the differen- 
tial cross section is reduced by about 30% (20%). For 
m(£+£-) > 130 GeV, the 0(a) QED corrections uni- 
formly reduce the differential cross section by about 12% 
in the electron case, and « 7% in the muon case. No 



4 Two Photon Radiation in W and Z Boson Pro- 
duction at the Tevatron 

The 0(a) EW corrections have a significant effect on the 
W and Z masses extracted by D0 and CDF. In particu- 
lar, the large shift in the masses caused by the emission 
of a photon from the final state lepton line raises the 
question of how strongly multiple photon radiation influ- 
ences the measured weak boson masses. At 0(a n ), W or 
Z decay with collinear emission of photons from a final 
state charged lepton gives rise to terms which are propor- 
tional to (a/7r)"ln"(A/2/ m 2) (y = W, Z) in n-photon 
exclusive rates. 

In order to find out how important multi photon ra- 
diation in W and Z production is for the measurement 
of Mw at the Tevatron, it is instructive to carry out a 
calculation of the two photon processes, pp — ► £1/77 and 
pp — » £ + £~ , yy. So far no calculation of these processes 
which is based on the full set of tree level 0(a 4 ) Feynman 
diagrams, and which is valid for arbitrary lepton-photon 
opening angles, has been carried out. For example, the 
calculation of Ref. Jl8| assumes that mt = 0. A non-zero 
ARi 7 cut, therefore, has to be imposed in order to avoid 
the collinear singularities. The Monte Carlo generator 
PHOTOS Ilia, on the other hand, treats final state pho- 
ton radiation in the leading-log approximation and thus 
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leads to results which can only be trusted in the collinear 
region. PHOTOS ignores initial state photon radiation 
altogether. 

In order to correctly take into account the effects of 
two photon radiation in W and Z production, a calcula- 
tion which gives correct results for small as well as large 
lcpton-photon opening angles is required. Here we report 
some preliminary results of such a calculation which is 
presently carried out E3. Due to the collinear singulari- 
ties associated with photon radiation from the charged 
lepton lines, there are many different peaks in the mul- 
tidifferential cross section. For an accurate evaluation 
of the cross section we therefore use a multiconfigura- 
tion Monte Carlo integration routine which automati- 
cally maps the peaks in the differential cross section to a 
uniform function according to the pole structure of the 
contributing Feynman diagramsE3. The matrix elements, 
taking into account finite kpton—masses, are calculated 
using the MADGRAPH package Eil, which automatically 
generates matrix elements in HELAS format£3. In order 
to maintain electromagnetic gauge invariance for £^77 
production in presence of finite W width effects, the W 
propagator and the WW-f and WW-f-f vertex functions 
in the amplitudes generated by MADGRAPH are modi- 
fied, using the prescription of Ref. ps| . 

In Table 0, we display the fraction of W — > ev and 
Z — ► e + e~ events (in percent) containing one or two 
photons at the Tevatron as a function of the minimum 
photon transverse energy, £^ in , for £^ in > 0.1 GeV, 
the approximate tower threshold of the electromagnetic 
calorimeters of CDF and D0. To simulate detector re- 
sponse, we have imposed the following acceptance cuts: 

pr(e) > 20 GeV, \q{e)\ < 2.5, and |?7( 7 )| < 3.6. (3) 

In the W case, we require in addition that 



j> T > 20 GeV 



(4) 



and 



65 GeV < M T (e + nj; v) < 100 GeV, (5) 

where Afi^e + Tvy; v) is the cluster transverse mass of the 
(e + wf)v system (n = 0, 1, 2). For Z events we require 

m(e + e-) > 20 GeV (6) 

and 

75 GeV < m(ee + 717) < 105 GeV, (7) 

where m(e + e~) (m(ee + rry)) is the e + e~ (ee + 717) in- 
variant mass. For muon final states, the fraction of 
events containing one (two) photons is roughly a fac- 
tor 2 (4) smaller than the results shown for W — > ev 
and Z — > e + e~. Table [l] demonstrates that a signifi- 
cant fraction of weak boson events contains two photons. 
Multiple photon bremsstrahlung thus is expected to have 
a non-negligible effect on the W mass extracted from ex- 
periment. 



Table 1: Fraction of W — > ev and Z — > e+e — events at the Tevatron 
(in percent) containing one or two photons with Ej-i^i) > E£} ln . 
Fractions are obtained by normalization with respect to the Born 
cross section. The cuts imposed are described in the text. 



E™ in (GeV) 


W — > evj 


W — » ei/77 


0.1 


28 


4.0 


0.3 


21 


2.3 


1 


14 


0.9 


3 


7.8 


0.2 


E^ in (GeV) 


Z — > e + e~7 


Z — > e + e~77 


0.1 


54 


16 


0.3 


42 


9.4 


1 


27 


3.9 


3 


15 


1.2 
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